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Abstract. For S a positive selfadjoint operator on a Hilbert space, 

d ^(t) + 2F(S)'^(t) + S 2 u(t) = 
dt dt 

describes a class of wave equations with strong friction or damping if F is a 
positive Borel function. Under suitable hypotheses, it is shown that 

u{t) = v(t) + w(t) 

where v satisfies 

2F(S)^(t) + S 2 v(t)=0 
dt 

and 

w(t) 

— > 0, as t — > +00. 

IK*)II 

The required initial condition v(0) is given in a canonical way in terms of u(0), 
u'(0). 



1. Introduction 

Let S be an injective nonnegative selfadjoint operator on a complex Hilbert 
space H. That is £ = S* > 0, ^ u p (S). Consider the damped wave equation 

(1.1) u"(t) + 2Bu'(t) + S 2 u(t) =0, t>0, 
with initial conditions 

(1.2) u(0) = f, u'(0) = g; 

here ' = d/dt. When B = 0, reduces to the wave equation, and the corre- 

sponding heat equation normally considered is 



■)'{t) + S 2 v(t) = 0. 
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In this paper we take B to be a positive selfadjoint operator which commutes with 
S and is " smaller than S" . More precisely, we assume that 

(1.3) = inf o-(S), 

i.e., is the spectrum of S = S* > 0, but is not an eigenvalue, F is a continuous 
function from (0, +oo) to (0,+oo), B — F(S), and F satisfies: there exists 7 > 
such that 

F(x) > x for < x < 7, 

F(rr) = 7 , 

(1.4) \ F(x) < x for x > 7, 
lim sup 2 ,_ >0 + F(x) < +00, 

liminf x _ >+00 ((l — 8)x — F{x)) > 0, for some S > 0. 

We make a further comment on (|1.3I) . Think of S* 2 as — A with suitable bound- 
ary conditions, acting on H = L 2 (fl) where f2 is a domain in R N . Then (jl.3|) implies 
that S 1 cannot have compact resolvent, and we are thus led to work exclusively in 
unbounded domains fl. 

The operator B represents a general friction coefficient. The most common 
case is the telegraph equation in which case 

B = al 

where a is a positive constant. In this case (|1.4j) holds with 7 = 0. Another simple 
case is 

B = aS a 

where the constants a, a satisfy 

a>0, ae[0,l). 

In this case, 7 = a 1 ^" in (|1.4j) . The only interesting case is when S is unbounded. 
The strongly damped wave equation refers to the case when B is also unbounded. 

Our main result, Theorem 13. 11 can be stated informally as follows. Let S,B = 
F(S) be as above and suppose /, g are such that (|1.1[) . (11.21) has a unique solution 
u. Consider the corresponding first order (in t) equation, obtained by erasing u"(t) 
in (|1.1[) and replacing u by v: 

(1.5) 2Bv'(t) + S 2 v(t) = 0, i>0, 
with initial condition 

(1.6) u(0) = h. 
This vector h is given by 

h = lx(oMS){( F (S) 2 - S 2 )HF(S)f + g) + /}, 

a formula which will be derived and explained later; and our conclusion requires 
that f,g are such that h ^ 0. We will show that 

(1 7) Ut) ■= ^- V ^ -> 
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as t — > 0, and we will find closed subspaces £ n of % such that £ n C £ n +i, 

oo 

I^J £ n is dense in H, 

71=1 

and 

l(t) < C„e- e "* 

for /, g G £ n where C n ,e n are positive constants. In general, there is no rate of 
convergence in (ll.2j) that works for all solutions. 

The point of the theorem is that, for large times the solution of the "hyperbolic 
equation" looks like the solution of the "parabolic equation" (|1.5|) . In the 

telegraph equation case when B = F{S) = al, (|1.5I) becomes 

2av'{t) + S 2 v{t) = 0, 

which, with (|1.6[) , is solved by 

v(t) = e~^ s2 h. 

In the case of strong damping, the solution of the limiting parabolic problem is 

v(t) = e~i B ~ ls2 h. 
Think of S = (-A)s on L 2 (R N ) and 

B = aS a =o(-A)f, 

< a < 1. Then B~ 1 S 2 = -(-A) 1 "^ with domain 

a 

D{B- 1 S 2 ) = H 2 - a (R N ), 

while 

L>(5 2 ) = ff 2 ^); 

here we use the standard Sobolev space notation. Thus _B _1 5 2 is a pseudodifferen- 
tial operator of lower order 2 — a than that of the Laplacian, unless a — in which 
case we have the telegraph equation. 

It has long been " known" that the telegraph equation is an asymptotic approx- 
imant for the heat equation, especially in the case of S 2 = —d 2 /dx 2 on L 2 (R). The 
pioneers in this area include G. I. Taylor |12j in 1922, S. Goldstein [§] in 1938 and 
M. Kac [9 in 1956. 

Some of the associated random walk ideas are discussed in [4] , which eventually 
led to pQ, the main theorem in which is the special case of our main theorem here 
with B — al. The importance of the interesting case of strong damping was 
recognized by Fritz Gesztesy and is discussed in detail in [5] . For additional results 
on strong damping, see [7]. 

Section [5] reviews some spectral theory. Section [3] contains the statement and 
the proof of our main result. The proof is patterned after that in pQ, but it in fact 
is simplified and streamlined. Section 2] contains examples. 
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2. Selfadjoint and normal operators 

Let S be a selfadjoint operator on H with spectrum a(S). By the spectral 
Theorem there exists an 1? space L 2 (A, S, v) and a unitary operator 

W : H -> L 2 (A,X» 

such that S is unitarily equivalent, via W , to the maximally defined operator of 
multiplication by a £— measurable function 

m : A -> er(S) C R, 

i.e., 

5/ = W" 1 M m W/, 

for 

/ £ £>(S) = {^-\g G H : mg e L 2 (A, S, i/)} 

and 

(M m5 f)(a;) = m(x) 5 (x), for x G A,g G L 2 (A,E,^). 
Two selfadjoint operators Si, S2 commute if and only if the bounded operators 
(A 1 /-5 1 )- 1 ,(A 2 /-5 2 )- 1 
commute for all Ai, A 2 6 C\l if and only if 

e ttSl and e lsS2 

commute for all t, s G M. Similarly, two normal operators Ni, N 2 with 

sup Re cr(Nj) < +00, j = 1,2 

commute if and only if e tNl and e sN2 commute for all t, s > 0; here N is normal 
means N — Si + iS 2 where Si, S2 are commuting selfadjoint operators. 

The functional calculus for S selfadjoint says that for every Borel measurable 
function F from a(S) CltoC, F(S) defined by 

F(S) = W- x M F(m) W 

is normal, and any two of these operators commute. Moreover, 

F -»• F(S) 

is linear and is an algebra homomorphism, thus 

Fi(S)F 2 (S) = (FiF 2 )(S), 

etc. Also, F(S) is bounded on H if and only if F is bounded on o~(S), and F(S) is 
selfadjoint if and only if F is real valued. And for S = S*, F(S) is semibounded 
(above or below) if and only if F(<r(S)) is, in M. 

In particular, for V a Borel set in [0, +00), Pr = Xr(S) is the orthogonal 
projection of H onto xr(S)(T-L); and 

P r F(S) = F(S)P r = P r F(SP r ) 

is the part of F(S) in T, and its spectrum is contained in V. 

If Fi, F 2 are complex Borel functions on a(S) that are bounded above, it follows 
that Fj(S) and T,^ =1 Fk(S) generate (Co) semigroups on % and 
(2.1) e «S=i*KS) =U n^ e tF k (S)^ 

and the product can be taken in any order. Finally, if L = F(S) = L* > 0, then 
[L]s denotes the unique nonnegative square root of L. 
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For more on the spectral theorem, see the books [6], [10] . [11] . 

3. The main result 

Consider the problem (|1.2[) . which we rewrite as 

(3.1) u" + 2Bu' + S 2 u = 0, t > 0, 

(3.2) u(0)=f,u'(0)=g, 
where S = S* > on H, 

(3.3) infer(S) = & a p (S), 

B = F(S) where F is a continuous function from (0, +oo) to (0, +oo) which is 
bounded near zero and strictly less than the identity function near infinity, in the 
sense that for some 5 > 0, 

(3.4) liminf((l - S)x - F(x)) > 0. 

x— v+oo 

We also assume there exists 7 > such that (|1.4j) holds, namely 
' F(x) > x for < x < 7, 



(3.5) 



F( 7 ) - 7, 

F(x) < x for x > 7, 

i limsup a .^ + F{x) < +00. 

Let r be the open interval (0,7) and let 

(3.6) P r = X r(S). 

Theorem 3.1. Assume all the statements in the above paragraph. Let v be the 
solution of the corresponding "parabolic" equation 

(3.7) 2Bv' + S 2 v = 0, 
obtained by deleting the second derivative term in (|3.1j) . Let 

(3.8) «(0) = h:= l -P T (f + [{B 2 - S 2 )P T ]^{Bf + g)). 

Then, for u the solution of (|3.ip . (|3.2I) . 

(3.9) u(t) = v(t)(l + o{t)) 

holds as t — > +00, provided h ^ 0. Moreover, if T n = [— , <S — — ] and i/0 ^ ft e 
Prni'H) for some n G N, £/ien 

(3.10) u(t) = u(t)(l +o(e- e "*)) 
/or some e n > 0. 

Remark 3.2. Note that 

Pr„(H) cP r „ +1 (H), 
an d U^°=i Pv n {T-L) is dense in Pt(H)- 
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PROOF. First recall that the square root of (B 2 — S 2 )Pr refers to the non- 
negative square root. We first show that the problem (|3.1|> . p.2[) is wellposed, by 
showing that it is governed by a (Co) contraction semigroup. 
We first treat the case of B = in (|3. II) . 

Svr 



Rewrite (13TT 



(1X21) as, for U = 
U' 



(3.11) 



Su' 





-s 



Su 



GU, 



Sf 
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U(0) = L = 

Let JC be the completion of D(S) © D(S) in the norm 

= (||m|| 2 + ||n|| 2 )^, 



where || • || is the norm on H. Then G defined by p. lip is skewadjoint on /C and 
thus generates a (Co) unitary group by Stone's Theorem (cf. [6], [10] . |llj). We 
consider this as a semigroup since we are only concerned with times t > 0. Next, 
(I3.1[) can be written in tC as 

U' = (G + P)U 

where 



P = 





-IB 



Since B = B* > 0, 



P 



||25n|| < (l-e)|[S , n[[+C e ||n| 



<(l-e) 



P 



for some e > and a corresponding C e > 0, for all n G D(S) C D(B), thanks 
to ([3TI)) and the last line in (153)) (or (TQl) ). Namely, write B = B x + B 2 := 



BP, 



(0,M) 



BP[m,oo)i where M is such that 



x > F(x) + Sx, 



i.e. 



F(x) < (1 - 5)x 

for x > M and F(x) is bounded in [0, M). Thus Bi is bounded, B 2 = B\ > and 

||B 2 n||<(l-*)||Sn|| + ||B 2 ||||n|| 



for all n G D(S). Thus for 



||P^V|k < (1-S)\\GN\\k + M\\N\\ k 

where 5 > and M = ||S 2 ||- It follows (cf. e.g. [6], [10], [11]) that G + P 
is rn— dissipative and generates a (Co) contraction semigroup on /C, since _P is 
obviously dissipative. Then (|3.1[) . p.2| has a unique strongly C 2 solution (resp. 
mild solution) if / G D(S 2 ),g G D(S) (resp. / G D(S),g G W): cf. [6] Chapter 2, 
Theorem 7.8]. 
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We shall express the unique solution using d'Alembert's formula. We seek 
solution of the form 



u(t) = e tc h 



where C is a Borcl function of S. By (|3.1I) . C must satisfy 

C 2 + 2BC + S 2 =0. 

Formally, 

C = C± = -B±(B 2 - 5 2 )i 
Selfadjoint operators have many square roots, but nonnegative selfadjoint operators 
have unique nonnegative square roots. Thus we uniquely define C± by 

(3.12) C± = -B ± (Q + %Q) 
where 

(3.13) Q = [(B 2 - S 2 ) X{ o,-r)(S)]KQ= l(S 2 - B 2 )\ Xh , +oo) (S)}? . 
Thus the solution u of (|3.1[) . (|3.2I) can be written as 

u(t) = e} c +h + + e tc -/i„, 

where C± are defined by (|3.12l) . (|3.13l) . There are strong C 2 solutions (resp. mild 
solutions) if and only if h± £ D(S 2 ) (resp. h± € D(S)). 

Given / = u(0),g — u'(0), we obtain h± by inverting the 2x2 system 

f = h+ + h_ 

g = C + h + + C-h_. 

An elementary calculation gives 

(3-14) h_ = ±(f-(Q + iQ)- 1 (Bf + g)) 

(3.15) h+ = \{f + (Qo + iQ)-\B.f + g)). 

Write 

where 



u = u\ + u 2 + u 3 

u 1 {t) = e tc +P (0 ^ ) h + , 
U 2 (t) = e tc +P [7 , +oo) /i+, 
u 3 (t) = e tc -/i_. 



First, 



IM*)II = \\e^ tQ e- tQo e~ tB h-\\ < \\e- tB h-\\ 
since e~~ lt ® is unitary and ||e~ t( 2 || < 1. Next, 

\\u 2 (t)\\ = || e ^e- tB P [7 , +oo) M < H^MI- 
The next estimate is the key one. 
For 

(3.16) h:=P (0 ^(h + ), 

1^)11 = \\e t( -- B+Qo) h+l 
We know that h € P(o )7 )(K): assume 

(3.17) 0^eP M (W)=:% 
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for some 8 > 0. Let 

Qos = QoP[8,j-S}- 

Then, since F(x) > x on [6, 7 — 5], F(x) — x > e on [6, j — 6} for some e > 0. Thus 

Qos > el. 

Consequently 

||t»i(t)||>e d ||e-* B fc||. 
It follows that for some constant Co, 

IMi)|| + ]M*)|| < C„e- rf |M*)ll- 

Thus 

(3.18) u{t) = Ul (t)(l + 0(e- et )). 

We must show that this holds with m replaced by v. 
The unique solution of (|3.7I) is 

v(t) = e-i B ~ ls \. 

Note that h as defined by (|3.8I) is Prh + where h + is as in (|3. 151) . To compare U\ 
with v, we need Taylor's formula with integral remainder, which for g G C 3 [0,/] 
and for some I > says that 

g(x) = ,g(0) + g'(Q)x + ^x 2 + I f " (x y) 2 g"'{y)dy. 



Applying this to 

(3.19) g(x) = 1- (1 -x)?, 0<x<l, 

yields 

(3-20) g ( L )f=l Lf + l I ?f + R f 

where R is a bounded operator commuting with L and satisfying 

R = R* > 0. 

Consequently 

\\ Ul {t) - v(t)\\ = || e -tfl(-/+fl- l Oo) ft _ e -* B ~ 13 *h\\ 

= \\ e -tB[I-(B- 2 Q 2 )i] h _ g-lB-^^U 

(3.21) = \\e-i B ~ ls2 {e-i L2 e- tR - I}h\\ 

by (USD, (O with L = B- 2 (B 2 - S 2 )P (0 , 7) = (I - B- 2 S- 2 )P (0n) . 
We have 

Cil < R < C2I 

on Us for some constants < Ci < C2 < +00. Furthermore, we also have 

CaJ < £ < 

on Hi for some positive constants £3, (4- It now follows from f|3 . 20[) that 
|| Ul (t) - v(t)\\ = \\e~i B ~ ls \l - e-i L e~ tR h)\\ 

and 

\\e-i L e- tR h\\ <e^\\h\\ 
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where 

C 5 = | + Ci>o. 

Consequently 

||tt(t)-«(t)||<||«(t)||0(e-^). 
Combining this inequality with p.l8j) yields the desired asymptotic relation 

\\u(t)-v(t)\\ _ t 
\\v(t)\\ ~ 0[e ] 

for some es > 0. 

Now let 7^ h € P( ,7)(^)- We must show that 

rq oo\ \Ht)-v(t)\\ A + 

(3.22) > 0, and i — > +oo. 

IK*) II 

We proceed by contradiction. Suppose p.22|) fails to hold for some h ^ in 
P(o, 7 )("H). Then, there exists ei > and — > +oo such that 

,„ O o, »"(*») ~^n)|| . 

for all neN. Choose 8 > and ft. G % = P[5. 7 _5]('H) (depending on ei) such that 

\\h-h\\ < j 

and let /, g be the corresponding initial data. Note that 

P[-y, + oo)l = P[-y, + oo)l 

for I = f,g, and / and g are modified only on the subspace -Pa(W) 

A := [5-5x,5 + 5i] U [7 - 8 - 81, 7 - 8 + Si], 

for some 81 > which can be chosen to be arbitrarily small. In particular, given 
£2 > we may choose /, g as above and additionally satisfying 

\\h\\ r 

M e[1 ~ £2 ' 1 + £2] - 

It follows that 

Kt)-«(t)||, [Ki)-5(t)||<^ 

for all t > 0. Consequently 



(3 24) IK*)II ||6(*)|| V:-e 2 y 4 

<™-«*(l±*) + e 4 61 



l-e 2 / 4 4 ' 

as t — !> +00, since ^ /i £ H<5, and To, £3 are positive constants depending on 8. 
But (|3.24[) contradicts f!3.23|) for < = £ n with n large enough. It follows that (|3.22p 
holds. This completes the proof of Theorem 13. II □ 



We remark that, in general, there does not exist a rate of convergence in (|3.22j) 
which works for all 5^ h E "H 7 . This follows from a careful examination of [l] and 
[2 Section 3]. 
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4. Examples 
Example 4.1. This is the simplest example. Take 
S 2 :=-A + w 2 I onL 2 {R N ) 
for w > 0. Similarly we can define 

S 2 k :=S 2k = (-A + w 2 I) k on L 2 (R N ), k G N. 
Theorem 13. II applies (with k G N, w > fixed) if we take 

B = aS% = aS ak 
for a>0, a£ [0,1). The resulting damped wave equation is 

(4.1) u u + 2a(-A + w 2 )^w t + (-A + w 2 ) k u = 0. 

This is a pseudo differential equation unless y e No = {0, 1, 2, ...}, in which case 
it is a partial differential equation. The simplest case of this is k = 3, a = |, in 
which case (|4.ip reduces to 

u u + 2a(-Au t + w 2 u t ) + (-A + w 2 ) 3 u = 0. 

As noted earlier, 7 is given by 

1 

7 = a 1 -" 

in the case of (|4.1I) , for all k G N. 

Example 4.2. We start with a brief summary of the "Wentzell Laplacian" 
discussed in [2j (cf. also [3]). Let fi be an unbounded domain in M. N with nonempty 
boundary dtt, such that for every R > there exists a ball B(xr,R) in D. Let 
A(a;) be an TV x N matrix for x € 51 such that A(x) is real, Hermitian and 

(4.2) ai|£| 2 <^K-£<«2|£| 2 
for all x e fi,£ G E w , where 

< ot\ < a-2 < 00 

are constants. Similarly let B(x) for x E dfl be a real Hermitian (iV — 1) x (N — 1) 
matrix satisfying (|4.2p with the same ai, a2- Assume A, B, and 7, /3 (introduced 
below) are sufficiently smooth. Define distributional differential operators on 51 
(resp. d(l) by 

Lu x = V • (A(ic)Vui), 
L a u 2 = V T • (B(x)V T u 2 ) 
for tii (resp. U2) defined on Jl (resp. dfl). Here V T is the tangential gradient on 
<90. The wave equation (without damping) we consider is 

(4.3) u t t — Lu in O, 

(4.4) Lu + (3d^u + 7 u - g/JLsu = on 89,. 
Here the conormal derivative term is 

d^u = {AVu) ■ v 

at x G 9f2, where v is the unit outer normal to dft at 1; £ C 1 (9il,R),/3 > 
0, 7 > 0, /3, ^ and 7 are bounded, and q G [0, +00). 

In [2], it is explained how the problem (14.31) . (|4.4[) can be rewritten as 

u" + S 2 u = 
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u(0) = /, u'(0)=g 

where the Hilbert space is 

n = L 2 {n,dx)®L 2 (dn,jj- 

Vectors in % are represented by u — {^f^j with u\ G L 2 (Q, dx) and u 2 G L 2 (dfl, . 
The norm in Ti. is given by 

h\\n = {IKIIi»(n) + \M\v>(aa,*r)}* ■ 

Here we write dT rather than the usual dS for the element of surface measure, since 
the letter S already is being used to denote the basic operator. 
The operator S 2 has the matrix representation 

-L 

^ y-qj3La y 

In [2] it was shown that S = [S 2 ] i , with a suitable domain, satisfies 
S = S* > 0, = inf a(S), g <r p (S). 

Furthermore, for all u G D(S), we have u = ( Wl ) , where U2 = tr(ui), the trace of 

_ ' _ W 
U\. As in Example l4.il Theorem 13. II applies to 

(4.5) u" + 2aS^u' + S 2k u = 0, 

u(0) = /, «'(0) = 9, fcGN. 

Again, this is a partial differential equation only when y€N. If k = 3 and a = |, 
the corresponding parabolic problem is 



S 2 = 



■>' + —S A v = 0, w(0) = h. 
2a 



The boundary conditions associated with (I4.5[) are 

Lw + fld^w + jw — qLgw = on <9f2 

for w = S 2 iu, j = 0, 1, fc — 1. 

Example 4.3 The simplest example of unidirectional waves in one dimension 
are described by the equation (for t, x G ffi) 

(4.6) itj = cua + k IIs: =: Mu, 

where (6, c) G R 2 \ {(0, 0)}. The most common case is c ^ 0, b — 0, in which case 
the corresponding equation for bidirectional waves is 

4 - M ) ( 4 + M ) u = u« - c 2 w ra = 0. 



v 9i / \dt 

The case of b ^ is the Airy equation, and (|4.6[) is the linearization of the KdV 
equation 

u t = cu x + bu xxx + ciuu x . 
For c = ^ i, the bidirectional version of (|4.6[) is 

M] ( + M ] u = u u — b 2 u xxxxxx = 0. 



0t \dt 
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Now, let % = L 2 (R), D = £ and T = -D 2 = T* > 0. Let 

S 2 = T 3 + a a T 2 + aiT, 
where ciq, a\ £ [0, +00). Consider 

0) = f(x), u t (x, 0) = g(x). 

In this case, 

B = aT = F(S) = F(T 3 + a T 2 + a x T). 
For x > 0, we want to consider the function 

G(x) = — (x 3 + a x 2 + a\x), 
a 

so that G(x) — x is negative in (0, 7) and positive on (7, 00) for some 7 > 0. But 

d ( G(x) — x\ d ( 1 . 9 . ,,\ 1. , 

— ) =-r( ~( x + a o* + (ai -a))}= -(2s + a > 0, 

aa;\a;/aa;\a /a 

and G(a;) = a; for a; 7^ if and only if 

x = ^(-00 ± \Z«o - 4(ai - a) ). 
Thus we get exactly one positive root if and only if 

(4.7) a>a 1 + ^, 

which we assume. It is now elementary to check that B — F(S) and F satisfies the 
assumptions of Theorem 13.11 In this case 

7 = l(~ a o + y^o-^ai-a) )■ 
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